Electron-positron pairs production in a macroscopic charged core 
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Classical and semi-classical energy states of relativistic electrons bounded by a massive and 
charged core with the charge- mass-radio Q/M and macroscopic radius Rc are discussed. We 
show that the energies of semi-classical (bound) states can be much smaller than the negative 
' electron mass-energy (— mc^), and energy-level crossing to negative energy continuum occurs. 

I Electron-positron pair production takes place by quantum tunneling, if these bound states 

are not occupied. Electrons fill into these bound states and positrons go to infinity. We 
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. explicitly calculate the rate of pair-production, and compare it with the rates of electron- 

CN ■ 

positron production by the Sauter-Euler-Heisenberg-Schwinger in a constant electric field. In 
addition, the pair-production rate for the electro-gravitational balance ratio Q/M = lO^i^ 
, is much larger than the pair-production rate due to the Hawking processes. 
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, ■ PACS numbers: 



INTRODUCTION 



uation for a relativistic electron 



•/^ ' As reviewed in the recent report \M , very soon after the Dirac ec 

^' WW n n 

Q ■ was discovered [2, l3[, Gordon [4] (for all Z < 137) and Darwin [5!] (for Z = 1) found its solution 

in the point-like Coulomb potential V{r) = —Za/r, they obtained the well-known Sommerfeld's 
formula for energy-spectrum, 
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I £{n,j) = mc? 



1 + 



Za 



(1) 



n-\K\ + (i^2 _ ^2^,2)1/2 

where the fine-structure constant a = /tic, the principle quantum number n = 1, 2, 3, • • • and 

^^|-(i + l/2) = -(/ + l), if J = ^ + i Z>0 
[ (j + 1/2) = 1, if j = / _ 1, / > 1 

/ = 0,1,2, - • • is the orbital angular momentum corresponding to the upper component of Dirac 
bi-spinor, j is the total angular momentum. The integer values n and j label bound states whose 
energies are £{n,j) £ (0, mc^). For the example, in the case of the lowest energy states, one has 



£{lSi) = mcVl - (Za)2, (3) 
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£(23,) = £{2P,) = mc^/l±v^l3W, (4) 
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S{2Pz_) = mc^sjl-^iZay. (5) 

For all states of the discrete spectrum, the binding energy mc? — £{n,j) increases as the nuclear 
charge Z increases. No regular solution with n = 1, / = 0, j = 1/2 and K = —1 (the lS'1/2 ground 
state) is found for Z > 137, this was first noticed by Gordon in his pioneer paper [Jj. This is the 
problem so-called "Z = 137 catastrophe" . 

The problem was solved [l-lj] by considering the fact that the nucleus is not point-like and 
has an extended charge distribution, and the potential V{r) is not divergent when r — )• 0. The 
Z = 137 catastrophe disappears and the energy-levels £{n,j) of the bound states 15, 2P and 
2S, • • • smoothly continue to drop toward the negative energy continuum [E- < —mc^), as Z 



increases to values larger than 137. The critical values Zcr for £{n,j) = —m(? were found 
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191 ]: Zcr — 173 is a critical value at which the lowest energy-level of the bound state 15*1/2 
encounters the negative energy continuum, while other bound states 2Pi/2; 25*3/2, • • • encounter the 
negative energy continuum at Zcr > 173, thus energy- level-crossings and productions of electron 
and positron pair takes place, provided these bound states are unoccupied. We refer the readers 



to 
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191] for mathematical and numerical details. 
The energetics of this phenomenon can be understood as follow. The energy-level of the bound 
state 15i/2 can be estimated as follow, 

iS(15i/2) = "T-c^ — < —mc^, (6) 

where f is the average radius of the 15i/2 state's orbit, and the binding energy of this state Ze^/f > 
2m(? . If this bound state is unoccupied, the bare nucleus gains a binding energy Ze^ /r larger than 
2mc^, and becomes unstable against the production of an electron-positron pair. Assuming this 
pair-production occur around the radius f, we have energies of electron (e_) and positron (e+): 

Ze" r-. — Ze" 



c|p_|)2 + m2c4-^; e+ = ./(c|p+|)2 + m2c4 + ^, (7) 

1^ V 'y 

where p± are electron and positron momenta, and p_ = — p+. The total energy required for a 
pair production is. 



e_+ = e_ + e+ = 2^(c|p_|)2 + m2c4, (8) 

which is independent of the potential V{r). The potential energies ^eV{f) of electron and positron 
cancel each other and do not contribute to the total energy ([8]) required for pair production. This 
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energy is acquired from the binding energy [Ze^ jf > 2m(P') by the electron filhng into the 
bound state l<S'i/2- A part of the binding energy becomes the kinetic energy of positron that goes 
out. This is analogous to the familiar case that a proton {Z = 1) catches an electron into the 
ground state 15'x/2) a-iid a photon is emitted with the energy not less than 13.6 eV. 

In this article, we study classical and semi-classical states of electrons, electron-positron pair 
production in an electric potential of macroscopic cores with charge Q = Z\e\, mass M and 
macroscopic radius Ra- 



il. CLASSICAL DESCRIPTION OF ELECTRONS IN POTENTIAL OF CORES 

A. Effective potentials for particle 's radial motion 

Setting the origin of spherical coordinates (r, 6, (p) at the center of such cores, we write the 
vectorial potential Af^ = {A,Aq), where A = and is the Coulomb potential. The motion 
of a relativistic electron with mass m and charge e is described by its radial momentum pr, total 
angular momenta and the Hamiltonian, 



H± = ±mc\ l + + (^)2 - Vir), (9) 

V mc mcr 

where the potential energy V{r) = bAq, and it corresponds for positive and negative energies. The 
states corresponding to negative energy solutions are fully occupied. The total angular momentum 
Ptf, is conserved, for the potential V{r) is spherically symmetric. For a given angular momentum 
Ptj) = mv±r, where v± is the transverse velocity, the effective potential energy for electron's radial 
motion is 



E±{r) = ±mc^/l + (^)2-y(r), (10) 
V mcr 

where ± indicates positive and negative effective energies, outside the core (r > Re), the Coulomb 
potential energy V{r) is given by 

Vondr) = — . (11) 
r 

Inside the core (r < Rc), the Coulomb potential energy is given by 

Ze^ 



Mn(r) 



(12) 



2Rc 

where we postulate the charged core has a uniform charge distribution with constant charge density 
p = Ze/Vc, and the core volume Vc = 47ri?|?/3. Coulomb potential energies outside the core (fTTj) 
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and inside the core (jl2p are continuous at r = Re- The electric field on the surface of the core, 

it J itc mc^Rc 

where the electron Compton wavelength Ag = h/{mc), the critical electric field Ec = m'^c'/{eK) 
and the parameter (3 is the electric potential-energy on the surface of the core in unit of the electron 
mass-energy. 



B. 



Stable classical orbits (states) outside the core. 



Given different values of total angular momenta p^, the stable circulating orbits Rl (states) are 
determined by the minimum of the effective potential Ej^[r) (|10|) (see Fig.[T|), at which dEj^[r) / dr = 
0. We obtain stable orbits locate at the radii Rl outside the core, 



Rl 



Ze^m j ^ 



1 



Ze2 

CPS 



, Rl > Rc 



(14) 



for different p^-values. Substituting Eq. into Eq. pOj) . we find the energy of electron at each 
stable orbit, 



£ = min(£'+) = mc^ 
For the condition Rl ^ Rc, have 



1 



Zf_' 
cp^^ 



(15) 



Ze^ 

cp^^ 



<1 
~ 2 



/3(4 + /32)l/2 _ ^: 



(16) 



where the semi-equality holds for the last stable orbits outside the core Rl ^ Rc + 0^ ■ In the 
point-like case Rc — ?• 0, the last stable orbits are 



cp^-^ Ze^ + i?L^O+, £-^0+. 



(17) 



Eq. (|15p shows that there are only positive or null energy solutions (states) in the case of a point- 
like charge, which corresponds to the energy-spectra Eqs. (13l4|5p in quantum mechanic scenario. 
While for » 1, radii of stable orbits Rl ^ ^ and energies £ — mc^ + 0~, classical electrons 
in these orbits are critically bound for their banding energy goes to zero. We conclude that the 
energies (llSh of stable orbits outside the core must be smaller than mc^, but larger than zero, 
£ > 0. Therefore, no energy-level crossing with the negative energy spectrum occurs. 
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C. Stable classical orbits inside the core. 

We turn to the stable orbits of electrons inside the core. Analogously, using Eqs. (1101120 and 
dE^{r)/dr = 0, we obtain the stable orbit radius Rl < 1 in the unit of Rc, obeying the following 
equation, 

P\Rl + .'Rl) = .^- n = ^^- (18) 
and corresponding to the minimal energy (binding energy) of these states 



Rr 



\ZeV Ri 2 



(3 - Rl) 



(19) 



There are 8 solutions to this polynomial equation (jlSp . only one is physical, the solution Rl that 
has to be real, positive and smaller than one. As example, the numerical solution to Eq. p8|) 
is Rl = 0.793701 for /3 = 4.4 • 10^^ and k = 2.2 ■ 10^^. In following, we respectively adopt 
non-relativistic and ultra-relativistic approximations to obtain analytical solutions. 

First considering the non-relativistic case for those stable orbit states whose kinetic energy term 
characterized by angular momentum term p^, see Eq. (llOh . is much smaller than the rest mass term 
mc^, we obtain the following approximate equation, 

^^Rl ^ ^\ (20) 

and the solutions for stable orbit radii are, 

^1/2 / r-T, \ 1/2 

and energies, 

£^(l-^P + ^/^Z?'/') mc\ (22) 
The consistent conditions for this solution are 

/3V2 > ^ for ji^ < 1, and /3 < 1 for non-relativistic 
limit v±_ <^ c, where the transverse velocity v± = p^/ {mRL). As a result, the binding energies (j22p 
of these states are m(? > > 0, are never less than zero. These in fact correspond to the stable 
states which have large radii closing to the radius Rc of cores and w j_ <^ c. 

Second considering the ultra-relativistic case for those stable orbit states whose the kinetic 
energy term characterized by angular momentum term p,^, see Eq. (llOp . is much larger than the 
rest mass term mc^, we obtain the following approximate equation, 

(3'Rl k\ (23) 
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and the solutions for stable orbit radii are, 



Rl 



1/3 



1/3 



(24) 



which gives Rl — 0.7937007 for the same values of parameters /3 and k in above. The consistent 
condition for this solution is /? > k ^ 1 for Rl < 1. The energy levels of these ultra-relativistic 
states are, 



Ze2 



1 



mc 



(25) 



and mc^ > £ > — 1.5/3mc^. The particular solutions £ = and £ ~ —mc^ are respectively given 



by 



Ze2 



1; 



(26) 



V^eV 3/?; 

These in fact correspond to the stable states which have small radii closing to the center of cores 
and v_\_ < c. 

To have the energy-level crossing to the negative energy continuum, we are interested in the 
values /3 > K ^ 1 for which the energy-levels (f25]) of stable orbit states are equal to or less than 



2/3 



1 



2 ^ 2 

mc < —mc . 



(27) 



\Ze^J 

As example, with /3 = 10 and k = 2, Rl ~ 0.585, <Sniin — — 9.87mc^. The lowest energy-level of 
electron state is •p^jiZc^') = n/ j3 ^ Q with the binding energy, 



f 

C-r, 



2^ 



(28) 



locating at Rl ~ {p^c/ Ze^y/'^ — )• 0, the bottom of the potential energy 14n(0) ([l2 



III. SEMI-CLASSICAL DESCRIPTION 



A. Bohr-Sommerfeld quantization 



In order to have further understanding, we consider the semi-classical scenario. Introducing the 
Planck constant fi = /i/(27r), we adopt the semi-classical Bohr-Sommerfeld quantization rule 



1, 



p^d(j)o=ih{l + -), 



p^{l)^n{l + -), / = 0,1,2,3,---, 



(29) 
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which are discrete values selected from continuous total angular momentum p^p in the classical 
scenario. The variation of total angular momentum Ap,^ = zizh in th unit of the Planck constant 
h, we make substitution 



21 + 1 



\Ze^J V 2Za 
in classical solutions that we obtained in section (IIIl). 



a 



(he) 



(30) 



1. The radii and energies of stable states outside the core (fHj) and become: 



Za 



2Za 



21 + 1 



£ = mc \ / 1 



2Za 
21 + 1 



(31) 
(32) 



where the electron Compton length A = h/{mc) 



2. The radii and energies of non-relativistic stable states inside the core ()2ip and (j22p become: 

'21 + 



Rl 
£ 



3. A(2Z + l)^i/2\ 2 



2Za 



(33) 
(34) 



3. The radii and energies of ultra-relativistic stable states inside the core (f24|l and (f25|l become: 

'2Z + l\i/3 



Rl 



2Za 



2/ + l\2/3 
2Za I 



(35) 
(36) 



Note that radii Rl in the second and third cases are in unit of Rc 



B. Stability of semi- classical states 



When these semi-classical states are not occupied as required by the Pauli Principle, the tran- 
sition from one state to another with different discrete values of total angular momentum / (/i,/2 
and A/ = ^2 — ^1 = il) undergoes by emission or absorption of a spin-1 [fi) photon. Following the 
energy and angular-momentum conservations, photon emitted or absorbed in the transition have 
angular momenta = p^pih) — P(i>{h) = ^{h — h) = and energy £^i = £{12) — £{h)- In this 
transition of stable states, the variation of radius is ARl = Riih) — RL{h)- 
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We first consider the stability of semi-classical states against such transition in the case of 
point-like charge, i.e., Eqs. (I31|32p with / = 0,1,2,- • •. As required by the Heisenberg inde- 
terminacy principle A(/>Ap(^ ~ 47rp^(/) ^ h, the absolute ground state for minimal energy and 
angular momentum is given by the / = state, ~ /i/2, Rl ~ \{Za)^^{l — i^Za)'^)^!'^ > 
and £ ~ mc^(l — {2Za)'^)^^'^ > 0, which corresponds to the last stable orbit (I17p in the classical 
scenario. Thus the stability of all semi-classical states / > is guaranteed by the Pauli principle. 
This is only case for Za < 1/2. While for Za > 1/2, there is not an absolute ground state in 
the semi-classical scenario. This can be understood by examining how the lowest energy states are 
selected by the quantization rule in the semi-classical scenario out of the last stable orbits (jlTh 
in the classical scenario. For the case of Za < 1/2, equating in Eq. (fT7j) to p^ = fi{l + 1/2) 
(j29p . we find the selected state / = is only possible solution so that the ground state / = in the 
semi-classical scenario corresponds to the last stable orbits (117[) in the classical scenario. While 
for the case of Za > 1/2, equating p^ in Eq. (fT7|) to p^ = h{l + 1/2) ([29]) . we find the selected 
semi-classical state 

i = ^>0, (37) 

in the semi-classical scenario corresponds to the last stable orbits (jl7p in the classical scenario. 
This state / = / > is not protected by the Heisenberg indeterminacy principle from quantum- 
mechanically decaying in ^-steps to the states with lower angular momenta and energies (corre- 
spondingly smaller radius (|3ip ) via photon emissions. This clearly shows that the "Z = 137- 
catastrophe" corresponds to i?^, — )• 0, falling to the center of the Coulomb potential and all semi- 
classical states (/) are unstable. 

Then we consider the stability of semi-classical states against such transition in the case of 
charged cores Rc ^ 0. Substituting p^ in Eq. ([29j) into Eq. (fT6|) . we obtain the selected semi- 
classical state / corresponding to the last classical stable orbit outside the core. 



Rc 
A 



' +1 -1 



-1/2 



p X 3/4 

(Za)i/4 ( :^ j > 0. (38) 



Analogously to Eq. (j37p . the same argument concludes the instability of this semi-classical state, 
which must quantum-mechanically decay to states with angular momentum / < I inside the core, 
provided these semi-classical states are not occupied. This conclusion is independent of Za-value. 

We go on to examine the stability of semi-classical states inside the core. In the non-relativistic 
case (1 » /3 > k^), the last classical stable orbits locate at i?L — and — )■ given by 
Eqs. (j21|22p . corresponding to the lowest semi-classical state ()33|34p with / = and energy mcP > 
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£ > 0. In the ultra-relativistic case (/? > k » 1), the last classical stable orbits locate at Rl — )• 
and — given by Eqs. (I24|25p . corresponding to the lowest semi-classical state ()35|36p with 
/ = and minimal energy, 

3 

mc^ ~ ——(^mc^. (39) 

This concludes that the / = semi-classical state inside the core is an absolute ground state in both 
non- and ultra-relativistic cases. The Pauli principle assures that all semi-classical states / > are 
stable, provided all these states accommodate electrons. The electrons can be either present inside 
the core or produced from the vacuum polarization, later will be discussed in details. 

We are particular interested in the ultra-relativistic case (3 > k ^ 1, i.e., Za ^ 1, the energy- 
levels of semi-classical states can be more profound than —mc^ {£ < —mc^), energy-level crossings 
and pair-productions occur if these states are unoccupied, as discussed in introductory section. 




IV. PRODUCTION OF ELECTRON-POSITRON PAIR 

When the energy- levels of semi-classical (bound) states £ < —mc^ (j27p . energy- level crossings 
between these energy-levels (j25p and negative energy continuum (jlOp for pr = 0, as shown in Fig.[2l 
The energy-level-crossing indicates that £ (j25l) and (llOp are equal. 



£ = E^, (40) 

where angular momenta in £ (j36p and (jlOp are the same for angular-momentum conservation. 
The production of electron-positron pairs must takes place, provided these semi-classical (bound) 
states are unoccupied. The phenomenon of pair production can be understood as a quantum- 
mechanical tunneling process of relativistic electrons. The energy-levels £ of semi-classical (bound) 
states are given by Eq. ([36]) or (1271). The probability amplitude for this process can be calculated 



by a semi-classical WKB method 



Wwkb(|p±|) = exp<j-|/^ Prdr\, (41) 



n 

where |p^| = p^/r is transverse momenta and the radial momentum. 



p^(r) = ^(c|p^|)2 + m2c4- + F(r)]2. (42) 

The energy potential V{r) is either given by Vout{r) ([TT]) for r > Rc, or Vin{r) (fT2|) for r < Rc. 
The limits of integration (j^T]) : Ri, = Ri < Rc ([2^ or ([55]) indicating the location of the classical 
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orbit (classical turning point) of semi-classical (bound) state; while another classical turning point 
Rn is determined by setting Pr{r) = in Eq. (|l2]l . There are two cases: i?„ < Rc and > Rc, 
depending on f3 and k values. 

To obtain a maximal WKB-probability amplitude (j4ip of pair production, we only consider the 
case that the charge core is bare and 

• the lowest energy-levels of semi-classical (bound) states: p^/[Ze^) = k//3 — t- 0, the location 
of classical orbit ([MD = i^fe ^ and energy <^ £ ^ £,nm = -3/3mcV2 (pHD : 

• another classical turning point i?„ < Rc, since the probability is exponentially suppressed 
by a large tunneling length A = Rn — Rb- 

In this case (i?„ < Rc), Eq. ()42p becomes 



and pr = leads to 

Rn / 2 \i/2 



Rc \l3m& 
Using Eqs. ()41l43l44p . we have 



:(c|p±|)' + mV]i/^ (44) 



W^WKB(IPJ) = expj ,^,(^,2^)1/2 X 

- exp i-0 8, 23/^[(c|pJ)^ + ^V]^/^i^e l 

Dividing this probability amplitude by the tunneling length A ~ Rn and time interval At ~ 
2TTh/(2nic^) in which the quantum tunneling occurs, and integrating over two spin states and the 
transverse phase-space 2 / dr±dp±/{2TTh)'^ , we approximately obtain the rate of pair-production 
per the unit of time and volume, 

d^N 1.15 /Za\ f 2.46 /iirX^/^' 



dtdSx ^ 67r2 Vri^sy^'^Pj (Za)i/2 A j j' ^^^^ 
1.15 f P \ r 2.46i?cl 



■ NS 



67r2 \tXR^J ^ I /3V2A 

where = Ze/Rc is the electric field on the surface of the core and the Compton time r = fi/m(? . 

To have the size of this pair-production rate, we consider a macroscopic core of mass M = Mq 
and radius Rc = 10km, and the electric field on the core surface Eg (fT3|) is about the critical field 
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{Es ~ Ec). In this case, Z = a-^{Rc/\f ~ 9.2 • lO^^, /? = Za\/Rc = Rc/X ^ 2.59 • W^^, and the 
rate (H7|) becomes 

dt(Px ~ 67r2 VrA3 7 V ^^^^ 1 "'^^ V A 



which is exponentially small for Rc S> A. In this case, the charge-mass radio Q/{G^/'^M) = 
2 • 10^^|e|/(Gi/^mp) = 8.46 • 10~^, where G is the Newton constant and proton's charge-mass radio 
\e\/{G^/^mp) = 1.1 • 10 
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It is interesting to compare this rate of electron-positron pair-production with the rate given 
by the Hawking effect. We take Rc = 2GM/c^ and the charge-mass radio Q/{G^/^M) ~ 
10~i^ for a naive balance between gravitational and electric forces. In this case /3 = 
\{Q/G^''^M){\e\/G^/'^m) ^ 10^, the rate (jUD becomes, 

Tns = S (^) (z^j) exp {-0.492(mM)} , (50) 



where the notation niM = Rc/ (2A). This is much larger than the rate of electron-positron emission 
by the Hawking effect 



Th ~ exp{-87r(mM)} , (51) 

since the exponential factor exp {— 0.492(mM)} is much larger than exp {— 87r(mM)}, where 
2mM = Rc/X > 1. 



V. SUMMARY AND REMARKS 



In this letter, analogously to the study in atomic physics with large atomic number Z, we study 
the classical and semi-classical (bound) states of electrons in the electric potential of a massive 
and charged core, which has a uniform charge distribution and macroscopic radius. We have found 
negative energy states of electrons inside the core, whose energies can be smaller than —mc^, and 
the appearance of energy-level crossing to the negative energy spectrum. As a result, quantum 
tunneling takes place, leading to electron-positron pairs production, electrons then occupy these 
semi-classical (bound) states and positrons are repelled to infinity. Assuming that massive charged 
cores are bare and non of these semi-classical (bound) states are occupied, we analytically obtain 
the maximal rate of electron-positron pair production in terms of the core radius, charge and 
mass. We find that this rate is much larger than the rate of electron-positron pair-production by 
the Hawking effect, even for very small charge-mass radio of the core given by the naive balance 
between gravitational and electric forces. 
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FIG. 1: In the case of point-like charge distribution, we plot the positive and negative effective potential 
energies E± ([TU|. p^/{mcRc) — 2 and Ze^ — 1.95mc^i?c, to illustrate the radial location (HH) of stable 
orbits where i?+ has a minimum ((T5|) . All stable orbits are described by cp^ > Ze^. The last stable orbits 
are given by cp^ — > Ze^ + 0+, whose radial location i?^ — > and energy £ ^ 0+. There is no any stable 
orbit with energy £ < and the energy-level crossing with the negative energy spectrum is impossible. 
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